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Abstract 


JWST is here. The early release observation (ERO) program provides us with the first look at the scientific data and 
the spectral capabilities. One of the targets from the ERO is HAT-P-18b, an inflated Saturn-mass planet with an 
equilibrium temperature of ~850 K. We present the NIRISS/SOSS transmission spectrum of HAT-P-18b from 0.6 
to 2.8 um and reveal the planet in the infrared beyond 1.6 um for the first time. From the spectrum, we see clear 
water and escaping helium tail features in an otherwise very hazy atmosphere. Our free chemistry retrievals with 
ATMO show moderate Bayesian evidence (3.79) supporting the presence of methane, but the spectrum does not 
display any clearly identifiable methane absorption features. The retrieved methane abundance is ~2 orders of 
magnitude lower than that of solar composition. The methane-depleted atmosphere strongly rejects simple 
equilibrium chemistry forward models with solar metallicity and a C/O ratio and disfavors high metallicity (100 
times) and a low C/O ratio (0.3). This calls for additional physical processes such as vertical mixing and 
photochemistry, which can remove methane from the atmosphere. 


Unified Astronomy Thesaurus concepts: Exoplanet atmospheres (487); Exoplanet atmospheric composition (2021); 


Exoplanet astronomy (486) 


Supporting material: machine-readable table 


1. Introduction 


The era of JWST has finally arrived (The JWST Transiting 
Exoplanet Community Early Release Science Team et al. 2022). 
The unparalleled infrared capability and photometric precision 
from JWST will evermore change our understanding of 
exoplanet atmospheres. With the expanded spectral windows 
into the infrared, we can now identify and measure new 
molecular features that have never been studied in detail before. 
Methane, an abundant chemical species on Earth and many solar 
system planets, has yet to be detected in a short-period exoplanet 
atmosphere. It is one of the most significant byproducts of life 
activities on Earth, and detecting methane is one of the many 
required steps to eventually finding evidence of carbon-based 
life as we know it beyond Earth. 

Methane is ample in solar system gas giants (Jupiter, Saturn, 
Uranus, and Neptune) where the upper atmosphere is cooler 
than 200 K. The lack of evidence for methane in short-period 
exoplanet atmospheres has been puzzling (Stevenson et al. 
2010; Kreidberg et al. 2019) as equilibrium chemistry forward 
models predict the rise of the methane abundance in sub- 
1000 K exoplanet atmospheres (Moses et al. 2011; Benneke 
et al. 2019). Methane can be destroyed at high temperatures 
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(Zahnle & Marley 2014) and it also has a short photochemical 
lifetime (Kasting et al. 1983; Miller-Ricci Kempton et al. 2012; 
Thompson et al. 2022). Measurements of methane abundances 
in exoplanet atmospheres will help us to better understand 
atmospheric vertical mixing based on interior thermal structure 
and photochemistry (Fortney et al. 2020). 

HAT-P-18b is a Saturn-like planet with a mass of 0.197 My, 
and a radius of 0.995 Rj,, (Hartman et al. 2012). It has an 
equilibrium temperature of ~850K. The relatively low 
temperature compared to other hot Jupiters (>1000 K) makes 
it an intriguing target for atmospheric characterization to study 
disequilibrium chemistry (Fortney et al. 2020), especially its 
impacts on methane abundance (Zahnle & Marley 2014; Tsai 
et al. 2021). 

Here we present the transmission spectrum of HAT-P-18b 
obtained using JWST Near Infrared Imager and Slitless 
Spectrograph single object slitless spectroscopy (NIRISS/ 
SOSS), which covers ~0.6 to 2.8 jum. Our data are sensitive to 
multiple molecular features from water and methane, which 
should allow us to robustly detect and measure their 
abundance. The long continuous wavelength coverage extend- 
ing into the optical also allows us to place tight constraints on 
HAT-P-18b's aerosol properties. In this Letter, we report the 
clear detection of water, an escaping helium tail, a sub- 
Rayleigh haze scattering slope, and methane depletion in the 
atmosphere of HAT-P-18b. 
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Figure 1. The top panel shows one of the integrations from rateints.fits files; the red lines show the first- and second-order spectra traces. The wavelength- 
dependent zodiacal light background and three zeroth-order nearby contamination stars are visible in the image. The bottom panel shows the reduced light curve at 
every pixel bin. Contamination 1 has a minimal effect on the planet spectrum as it is far away from the first order and not extracted in the second order. We detect a 
time-varying flux level in contamination 2 (~12% dimming), which is likely caused by an eclipsing binary. Contamination 3 has a spike going across the second-order 
spectrum, which leads to wavelength-dependent dilution effects. We removed the contaminated parts of the spectrum from sources two and three in the analysis. 


2. Observations and Data Reduction 


This data set was collected as part of the early release 
observation program (program ID: 2734; Pontoppidan et al. 
2022). A transit of HAT-P-18b was observed using NIRISS/ 
SOSS on 2022 June 13 with the GR700XD grism and the 
CLEAR filter. The SUBSTRIP256 detector was used to capture 
both order one and order two of the dispersed spectra. A total 
of 7.15 hr of time series observation was taken to ensure 
sufficient out-of-transit baseline coverage. There are nine 
groups per integration with each integration taking 54.94 s 
and a total of 469 integrations. 


2.1. Data Reduction 


We begin our data reduction with the rateints.fits 
data downloaded from the MAST archive on 2022 July 28. The 
products are generated from the JWST Science Calibration 
Pipeline stage 2a, which fits the accumulating signal ramp for 
every pixel in each integration. One integration of the 
rateints.fits file is shown in the top panel of Figure 1. 
The image from the SUBSTRIP256 subarray has a dimension 
of 2048 x 256 pixels with a wavelength increasing from right to 
left. Since SOSS mode is slitless, contamination from nearby 
stars and the sky background can dilute transit signals. The first 
step is to remove the zodiacal light background, which increases 
with wavelength with a sharp drop around 2.1 jum caused by the 
edges of the pick-off mirror. We used the smoothed background 
template from the commissioning observation (observation 5 of 
program 1541) downloaded from jwst-docs. The back- 
ground template is first scaled based on the median of the 
bottom 20 percentile pixels between columns 300 and 500 
(72.3 to 2.5 um) and then subtracted from every integration. 
Since we do not expect the zodiacal light background to be 
time-dependent on our 7.15 hr timescale, we subtract the same 
scaled background template from all integrations. 


2 https: / /jwst-docs.stsci.edu/ 


Next, we determine the first- and second-order spectral trace 
locations by first taking the point-spread function from the bright 
first-order spectral region (~1.5 jum) and then cross correlating it 
at every column. The peak cross-correlation values at each 
column as a function of column pixel values are then fitted with a 
polynomial to smooth out the effects of contamination and bright 
outlier pixels. The first-order spectrum is then extracted with a 
width of 60 pixels centered at the trace as shown in red lines in 
the upper panel of Figure 1. The scipy.interpolate. 
interpld function is used at each column to interpolate each 
extraction at the subpixel levels. The same steps are then used to 
extract the second-order spectra after the first order is extracted 
from the 2D image. 

After the first- and second-order spectra have been extracted 
and aligned in the y-axis, we applied a median filter along each 
row to remove all outlier bright pixels within the 60 pixel wide 
spectral cut. Then we subtract the median of the combined top 
and bottom 15 pixels at each 60 pixel wide column to remove 
the 1/f noise. At last, we sum the middle 30 pixels at each 
extracted spectral column for every integration to get the light 
curves for all wavelength channels. All the cleaned light curves 
at each column are shown in the lower panel of Figure 1. 


2.2. Light-curve Fitting 


We first fit the white light data using the Batman 
(Kreidberg 2015) transiting light-curve model with emcee 
(Foreman-Mackey et al. 2013). There are a total of six free 
parameters in the fit including the midtransit time, transit depth, 
semimajor axis, inclination, linear slope, and a constant offset. 
There is a starspot crossing event during the transit that changes 
the transit light-curve shape in a wavelength-dependent way 
(left panel of Figure 2). We mask out the starspot crossing 
Gntegration 240 to 272) during the white light fit. After 
subtracting the best-fit model, we fit a polynomial model to the 
starspot crossing residual. Since the starspot shape varies as a 
function of wavelength due to its different temperature relative 
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Figure 2. The transit light curve has a clear starspot crossing event starting around midtransit (left panel). The starspot crossing affects the light curve as a function of 
wavelength with shorter wavelengths being the most affected. We first measure the amplitude of the starspot at each wavelength and then use the PHOENIX stellar 
models to obtain a best-fit starspot temperature of 3596 + 31 K (right panel). Assuming other unocculted starspots are at a similar temperature, we then use the stellar 
flux variability level seen in TESS observations to estimate the dilution effect. We obtain a «20 ppm level dilution effect between 0.6 and 2.8 um, which is much 


smaller than the uncertainties in the transmission spectrum. 


to the stellar photosphere, we use the best-fit model from the 
white light fits and scale it for each wavelength (Sing et al. 
2011). 

The best-fit white light midtransit time (2459743.853395 + 
0.000023 MJD), semimajor axis (a/Retar = 16.52 + 0.06), and 
inclination (88.66 + 0.03 degrees) are then fixed to fit the light 
curve from each column. There are four free parameters in the 
fit including the transit depth, a linear slope, a constant offset, 
and a scaling factor for the starspot crossing model. We used 
the scipy.optimize.curvefit function for individual 
wavelength channel fits. For the limb darkening, we used 
quadratic coefficients from the 3D stellar model in the Stagger 
grid (Magic et al. 2015). For the wavelength solution, we used 
the jwst-niriss-wavemap-0021.fits wavelength map 
file from the JWST calibration reference data system. 

We compared our final reduced transmission spectrum to an 
independent reduction from Néstor Espinoza using the pipeline 
TransitSpectroscopy” and they are in excellent agree- 
ment as shown in the top panel of Figure 4. 


2.3. Contamination and Dilution 


To ensure the accuracy of our transmission spectrum we 
have investigated possible sources of contamination and 
dilution. We identify the following three sources that could 
potentially affect the planet spectrum. 


2.3.1. Zeroth-order Contamination 


As NIRISS/SOSS is a slitless spectrograph, nearby stars can 
overlap with the target’s spectral trace. There are three clear 
zeroth-order nearby star contaminations in the observation as 
shown in the upper panel of Figure 1. Contamination 1 is on 
the second-order spectral trace but fortunately, it only affects 
the longer wavelength portion of the second order, which is not 
extracted due to the relatively low counts. It is also too far away 
from the first-order spectrum to cause dilution. Contamination 
2 is more complex as it sits directly on the first-order dispersion 
and its brightness also varies in time (Figure 1). It is likely an 


a https:/ /github.com/nespinoza/transitspectroscopy 


eclipsing binary with ingress starting ~2.4 hr after HAT-P- 
18b’s midtransit; the observation only catches a partial egress 
but it has a V-shaped transit with a maximum dimming of 
~12%. We masked this part of the spectrum (1.98-2.01 um) in 
the analysis. Contamination 3 is located between the first and 
second order, but it has a long tail going across the second 
order around 0.7 um. The tail leads to a wavelength-dependent 
dilution in the transit depth; therefore we removed the 0.68 to 
0.72 um part of the spectrum in the analysis. 


2.3.2. A Nearby Source from Gaia 


There is a nearby star (2/66 away) in the Gaia database (Gaia 
Collaboration et al. 2018; source ID: 1334573817793501696) 
that could cause dilution in the overlapped spectral dispersion. 
The companion has a delta magnitude of 4.587 and 4.898 fainter 
than HAT-P-18b, respectively, in the Gaia blue and red 
passbands. We performed a Gaia colors fit of the companion 
relative to HAT-P-18 (Teg = 4800 K ) and obtained a best-fit Tore 
of 5500 K + 300 K. Then we normalize the 5500 K and 4800 K 
PHOENIX models to the delta magnitude and calculate the 
dilution effect on the transit depth. Assuming the worst-case 
scenario of complete companion spectral overlap, the resulting 
dilution is less than 20 ppm between the 0.6 and 3 um regions of 
the transit spectrum. This is much smaller than the uncertainties 
on the final transmission spectrum, and we only expect slight 
spectral overlap from the companion base depending on the 
position angle of the instrument. Therefore, we conclude it is safe 
to ignore this Gaia nearby source. We also checked archival 
infrared high-contrast imaging data and there are no other known 
nearby sources (Ngo et al. 2015, 2016). 


2.3.3. Unocculted Starspots 


Unocculted starspots can lead to a wavelength-dependent 
dilution effect on the transmission spectrum (Agol et al. 2010; 
Pont et al. 2013). To estimate the magnitude of this effect, we 
follow the steps detailed in Section 3.5 of Sing et al. (2011). 
We first fit for the starspot temperature using the best-fit 
starspot scaling factor at each wavelength. The best-fit starspot 
temperature is 3596 + 31 K assuming a host star temperature of 
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4800 K (Figure 2). Next, we estimate the photometric stellar 
variability amplitude from TESS observations, which observed 
the star during sectors 25 and 26. After removing the transits in 
TESS, we used a Lomb-Scargle periodogram to find a best-fit 
variability period of ~14.3 days, which is consistent with the 
reported stellar rotational period of 14.66 + 0.03 days (Everett 
et al. 2013). The corresponding variability amplitude is 
0.072496 and it is within the 0.1% level mean R-band 
brightness variation observed from the 12 month ground-based 
monitoring campaign (Kirk et al. 2017). Finally, we use 
Equations (4) and (5) from Sing et al. (2011) to calculate an 
expected dilution effect of «20 ppm between 0.6 and 3 um in 
the transmission spectrum. Since this effect is much smaller 
than the uncertainties on the transmission spectrum, we 
conclude it is safe to ignore the unocculated starspots. 


3. Detection of Excess Helium Absorption and Tail 


We detect excess metastable 1083 nm helium absorption in the 
transmission spectrum. The helium line is not fully spectrally 
resolved with the SOSS (R ~ 550 at 1083 nm) but the transit 
depth from a 2 pixel wide bin (one resolution element) centered 
at 1083.065 and 1083.998 nm is visibly deeper than that in 
nearby pixels (Figure 3). We measure a helium transit depth of 
1.97% + 0.016% with the 2 pixel bin relative to a nearby white 
light (40 pixel bin) depth of 1.88% + 0.004%. This is consistent 
within 1c to the 2.46% + 0.15% helium transit depth reported in 
Vissapragada et al. (2022) and Paragas et al. (2021) after 
accounting for the spectral resolution difference. Although our 
measurement is more diluted than that of Vissapragada et al. 
(2022), our detection of excess absorption is more statistically 
significant (5.60 versus 4.30). 

We also see evidence of a helium tail at egress (Figure 3). 
We subtracted the nearby white light from the helium light 
curve and the residual clearly shows asymmetry between 
ingress and egress. Future detailed follow-up 3D hydrodyna- 
mical simulations of the tail may tell us more about the 
escaping material and stellar environment of HAT-P-18b. A 
planetary tail from escaping helium has only previously been 
detected on WASP-107b (Allart et al. 2019; Kirk et al. 2020; 
Spake et al. 2021), WASP-69b (Nortmann et al. 2018), and 
HD 189733b (Guilluy et al. 2020). Previous ground-based 
observations of HAT-P-18b did not see the tail because the 
signal-to-noise of their excess absorption detection was lower 
than that of the detection we report here. This demonstrates the 
unmatched photometric precision from JWST compared to 
typical ground-based observatories. Although lower in spectral 
resolution compared to ground-based high-resolution spectro- 
graphs (R > 10000), JWST NIRSpec/G140H (R — 2700) is a 
promising instrument mode for helium observations thanks to 
its superior stability and precision. 


3.1. Escaping Helium Modeling 


We fit the metastable helium (Seager & Sasselov 2000) 
transmission spectrum between 1.07 and 1.10 jum with a 1D, 
isothermal Parker-wind model (Parker 1958) based on the p- 
winds code (version 1.3.0; Dos Santos et al. 2022) to 
determine the mass-loss rate of HAT-P-18b. This code applies 
the original formulation from Oklopéié & Hirata (2018) along 
with modeling improvements proposed in Lampón et al. (2021) 
and Vissapragada et al. (2022). We assume the planetary 
outflow has a He/H fraction of 10/90. We further assume that 
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Figure 3. We detect escaping helium in HAT-P-18b. The top panel shows the 
light curve centered at the 1083 nm helium line vs. white light light curve from 
nearby pixels. The middle panel shows the residual (binned in red) from 
subtracting the outside helium light curve from the inside helium light curve. 
The residual clearly shows a deeper transit and a tail after the egress. The 
vertical gray dashed lines show the first and fourth contact and the horizontal 
gray dashed line is the zero residual level. The bottom panel is the transmission 
spectrum of HAT-P-18b at pixel level overplotted with the best-fit p-winds 
model. 


the high-energy spectral energy distribution of HAT-P-18 is 
similar to HD 85512, and adopt the MUSCLES spectrum 
(Loyd et al. 2016; Youngblood et al. 2016) of the latter as a 
proxy when calculating the ionization balance in the planetary 
outflow. We took into account the Roche lobe effects and the 
fact that the planet is only irradiated over 7 steradians rather 
than the 47 steradians assumed by p-winds models. 

Our forward models indicate that the mass-loss rate of 
3.7 x 10? g s^! estimated by Vissapragada et al. (2022) for an 
outflow temperature of 8000 K is consistent with the observed 
JWST transmission spectrum in the metastable He feature (see 
Figure 3). Since the absorption profile is not resolved by 
NIRISS, there is a degeneracy between the escape rate and 
outflow temperature when fitting isothermal Parker-wind 
models to the data. Because of this degeneracy, solutions with 
lower mass-loss rates will fit the observation just as well if they 
have lower temperatures. This effect is also seen in the narrow- 
band photometry technique used by Vissapragada et al. (2022), 
and it is likely going to be a limitation in future NIRISS 
observations of He in other planets as well. High-resolution 
spectroscopy from the ground can resolve the planetary 
absorption and break this degeneracy. 


THE ASTROPHYSICAL JOURNAL LETTERS, 940:L35 (8pp), 2022 December 1 


0.0210 


0.0205 


0.0200 


0.0195 


0.0190 


Transit Depth 


0.0185 


0.0180 


0.0175 
0.0210 


0.0205 


0.0200 


0.0195 


0.0190 


Transit Depth 


0.0185 


0.0180 


0.0175 


0.002 


0.001 


0.000 


Residual 


—0.001 


—0.002 


è Fu 
* Espinoza 
Kirk (ground)+ Tsiaras (WFC3) 


PHOENIX forward model (Z=Zo C/0=0.53) 
PHOENIX forward model (Z=100xZ_ C/0-0.3) 


ATMO retrieval free chemistry without CH4 
ATMO retrieval free chemistry 


1.0 LS 


BIC=1772.39 
BIC-1199.06 


BIC-1187.21 
BIC-1190.19 


2.0 


Wavelength (um) 


Xi1.715 
x2=1.149 


log evidence=5570.77 
log evidence=5574.56 


2.5 


Fu et al. 


Figure 4. The top panel shows the comparison between two independent data reductions and previous observations from the ground (Kirk et al. 2017) and HST/ 
WFC3 (Tsiaras et al. 2018) with no offsets applied. All data show excellent consistency with each other. The middle panel shows the spectrum (Fu version) 
overplotted with two PHOENIX forward models and two ATMO free chemistry retrieval best-fit models. The higher R — 700 spectrum is in gray and the lower 
R — 200 spectrum is in black. The retrieval is performed on the R — 700 spectrum. The lack of a clearly identifiable methane feature rules out the solar composition 
atmosphere. The 100 times solar metallicity and 0.3 C/O forward model provide a decent fit but underpredict the 1.4 um water feature while overpredicting the 
2.3 uum methane feature. The ATMO free chemistry retrieval (with H20, CO2, CO, and CH4) provides the best-fit model, and an identical retrieval without CH, gives a 
comparable fit. The log Bayesian evidence difference of 3.79 corresponds to moderate evidence supporting the detection of methane. The BIC metric slightly favors 
the retrieval without methane due to one less free fitting parameter. The bottom panel shows the corresponding residuals from the models in the middle panel. All 


models detect clear water absorption and a hazy atmosphere. The escaping helium part of the spectrum was excluded and modeled separately. 


4. Atmospheric Retrievals and Modeling 
4.1. PHOENIX Forward Model 


PHOENIX is a self-consistent atmospheric forward model 
assuming local thermodynamic equilibrium and equilibrium 


chemistry (Lothringer et al. 2018). The transmission spectrum 
shows distinct haze (optical scattering slope), water, and the 


lack of strong methane features. Solar composition PHOENIX 


models fail to fit the spectrum as they largely overpredict the 
methane features (Figure 4). One way to suppress the methane 
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Figure 5. We explored the effects of vertical mixing and photochemistry using VULCAN (Tsai et al. 2021). The blue line is chemical abundance for H20, CH4, CO, 
and CO, at each pressure level from the PHOENIX forward model assuming equilibrium chemistry and solar composition. The orange, green, and red lines show the 
abundance with different K,, values. The dashed lines are with photochemistry turned on. The dashed black line and gray shaded regions are the retrieved constraints 
from ATMO. We see both vertical mixing and photochemistry are effective at removing methane from the atmosphere. 


while preserving water is to increase the metallicity and 
decrease the C/O ratio. Therefore we generate a PHOENIX 
forward model grid at 10, 50, and 100 times solar metallicity 
and C/O ratios of 0.3 and 0.53 with various haze strengths and 
slopes. The best-fit PHOENIX model oe = 1.149 with 1013 
degrees of freedom) has a 100 times solar metallicity and a 0.3 
C/O ratio as shown (red) in the midpanel of Figure 4. It 
underpredicts the 1.4 jum water feature while overpredicting 
the 2.3 um methane feature. This suggests that the equilibrium 
chemistry model struggles to effectively remove methane 
without also oversuppressing the water features. 

We further explored methane depletion by running models 
with molecular abundances quenched at 0.1, 1 bar to simulate the 
effects of vertical mixing. The methane abundance at 0.1 bar is a 
few orders of magnitude below its maximum abundance, 
resulting in the most methane depletion if quenching occurs 
near this level in the atmosphere (Figure 5). We then pick the 
0.1 bar as the quench point and were able to reduce the best-fit 
forward model metallicity from 100 to 50 times solar. The 
0.1 bar pressure level was selected solely based on examining the 
methane abundance profile and finding the corresponding 
pressure level for the minimum methane mixing ratio. Only 
CH4, CO;, CO, and H5O were quenched in our model and 
quenching approximation does not work for every chemical 
species (Tsai et al. 2017). 

We also tested the effects of photochemistry and vertical 
mixing using the chemical kinetics code VULCAN (Tsai 
et al. 2021). We ran VULCAN using the atmosphere structure 
from the best-fit PHOENIX model with default parameters 
using the HCNO chemical network with various vertically- 
constant K,, both with and without photochemistry. Abun- 
dance profiles of H20, CH4, CO, and CO, are shown in 
Figure 5. Various scenarios with vertical mixing and 
photochemistry can lead to significant methane depletion, 
matching constraints from the retrieval, but determining the 
specific contribution from each effect at work in HAT-P-18b 
requires a more detailed retrieval analysis, which is beyond of 
the scope of this Letter. 


4.2. ATMO Free Chemistry Retrieval 


To better explore the abundance and significance of the 
detected chemical species, we also ran a free chemistry retrieval 
with ATMO (Tremblin et al. 2015, 2016) including temper- 
ature, radius, cloud deck pressure, haze strength, haze 
scattering slope, H20, CO», CO, and CHy. We detected water 
with a volume mixing ratio (VMR) of log(H5O) = —3.03 031) 


and retrieved a methane abundance of log(CH4) = —5.08703;. 


Compared to the PHOENIX forward model, the water 
abundance is consistent with a ~5 times solar metallicity 
atmosphere at 1 mbar assuming chemical equilibrium and a 
solar C/O ratio while methane is about 2 orders of magnitude 
lower. We find a CO abundance of log(CO = —4.76*15$3) and a 
CO, abundance of log(CO; = —4.19704) are more weakly 
constrained due to the lack of spectral features at these 
wavelengths, but their low retrieved abundances are consistent 
with the forward model predictions since most of the carbon 
should be in the form of methane. 

To further investigate the presence of methane in the 
atmosphere we ran another identical retrieval but without CH4. 
The best-fit retrieval model (Bayesian information criterion, 
BIC = 1190.19 and log Bayesian evidence = 5570.77) pro- 
duced a comparable fit to the first full retrieval (BIC=1187.21 
and log Bayesian evidence — 5574.56) with methane included. 
The difference of 3.79 in log Bayesian evidence would 
correspond to moderate evidence (Trotta 2008) supporting 
the detection of methane. However, the BIC favors the retrieval 
without methane due to one less free fitting parameter. The 
retrieved VMR of water (H20) = —3.8*('/3) stays unchanged 
and CO, (log(CO2) = 2572758 remains unconstrained, but 
the CO (log(CO) = —3.387291) abundance increases. We 
believe this is due to the overlapping feature from CH, and 
CO around 2.3 jim. Based on the free chemistry retrieval 
results, we are hesitant to conclude clear detection of methane 
in the atmosphere of HAT-P-18b. A definitive detection and 
more stringent limits on CH4, CO, and CO; could be made with 
NIRSpec G395H, as the molecular signatures of all three 
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Table 1 
HAT-P-18b Transmission Spectrum 

Wavelength Midpoint Transit Depth Uncertainty 
(um) (ppm) (ppm) 
0.6078 18908 1009 
0.6087 18903 616 
0.6096 19277 598 
0.6104 19544 653 
0.6113 19468 525 


(This table is available in its entirety in machine-readable form.) 


carbon-bearing species are considerably stronger at longer 
wavelengths. 


5. Discussion and Conclusion 


The lack of identifiable methane features in the transmission 
spectrum (Table 1) challenges the simple equilibrium chemistry 
assumption and motivates the exploration of additional 
physical processes such as vertical mixing and photochemistry 
(Cooper & Showman 2006; Steinrueck et al. 2019; Drummond 
et al. 2020). Vertical mixing brings up material from deeper in 
the planet and changes the atmospheric composition probed by 
the transmission spectrum (Zahnle & Marley 2014; Fortney 
et al. 2020). Depending on the planet's thermal structure and 
horizontal circulation, mixing could lead to the destruction of 
methane from the hotter interior or dayside of the planet 
(Komacek et al. 2019). On the other hand, photochemistry at 
the top of the atmosphere can break down methane into other 
hydrocarbon and haze precursors (Tsai et al. 2021). The 
detection of a strong haze layer on HAT-P-18b is consistent 
with the presence of a methane-derived photochemical haze. 

We present the NIRISS/SOSS transmission spectrum of 
HAT-P-18b, an inflated Saturn-like mass planet with an 
equilibrium temperature of ~850K. The 0.6 to 2.8 um 
wavelength coverage from SOSS gives us the first look at this 
planet in the infrared. We detect water, an escaping helium tail, 
and methane depletion in the hazy atmosphere of HAT-P-18b. 
The water and helium spectral features are consistent with 
previous observations (Tsiaras et al. 2018; Paragas et al. 2021; 
Vissapragada et al. 2022). The detection of a helium tail 
showcases the unparalleled level of photometric precision from 
JWST compared to ground-based observation. JWST NIR- 
Spec/G140H (R ~ 2700) will be the ideal instrument mode for 
any future study of escaped helium. 

The methane-depleted atmosphere adds to the growing 
observational evidence (Stevenson et al. 2010; Benneke et al. 
2019; Kreidberg et al. 2019) of missing methane in sub-1000 K 
planets. This calls for including disequilibrium chemistry 
processes such as vertical mixing and photochemistry, which 
are effective at delaying the onset of detectable methane features 
to lower equilibrium temperatures. The fast-growing JWST data 
sets of exoplanet atmospheres will certainly challenge our simple 
model assumptions and additional modeling efforts with more 
physical processes will be needed. 

Future follow-up observation of HAT-P-18b with JWST 
NIRSpec/prism covering the 2.8 to 5.3 um wavelength range 
will be useful at determining how depleted methane is in the 
atmosphere by measuring the stronger 3 to 4 jum methane 
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feature. Additionally, it will better constrain the planet’s C/O 
ratio by simultaneous coverage of CO and CO, features at 4 
and 5 um. 
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